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Abstract 

We present a computation of QCD next-to-leading order virtual corrections to the 
top production and decay process at linear colliders. The top quarks are allowed to be 
off-shell and the production and decay subprocesses are treated together, thus allowing 
for interference effects. The framework employed for our computation is the double 
pole approximation (DPA). We describe the implementation of this approximation 
for the top production and decay process and compare it with the implementation 
of DPA for the evaluation of QED corrections to the W pair production at LEP II. 
Similarities and differences between the two cases are pointed out. One result of interest 
is the incomplete cancelation of interference corrections. Other results include values 
for the total NLO top production cross section, and the impact the nonfactorizable 
(interference) corrections have on the top invariant mass distribution. 

1 Introduction 

The top quark might prove to be one of the most interesting elementary particles discovered 
so far. Since its mass is so large (close to the electroweak symmetry breaking scale), it can 
be assumed that study of top properties will allow us to learn new things about the physics 
at the next energy scale. So far, though, due to low statistics in the Run I at Tevatron, our 
knowledge of the top quark is rather limited. Future runs at the Tevatron and the LHC will 
provide more information, but the analysis of this new data will be complicated by large 
QCD backgrounds. 

An e^e“ collider with center-of-mass energy at and above the top threshold promises 
to provide a clean environment in which to perform precision studies of the top quark [1]. 
Indeed, it is conceivable that at such a machine the study of top can be performed with a 
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precision similar to that aquired in the study the W boson at LEP II. This means order % 
(and better) measurements of differential cross sections for processes involving the top quark. 
Such a precision in measuring experimental quantities implies the need for a like precision 
in our theoretical understanding of these processes. This in turn requires the inclusion of 
radiative corrections in our predictions. 

A wealth of information about the top quark (like its mass, width, strong and Yukawa cou¬ 
pling constants) can be extracted from measurements performed at the production threshold. 
For this region, comprehensive theoretical studies (NNLO computation with resummation 
of large logarithms, careful treatment of the renormalon ambiguity) have already been per¬ 
formed ([2] and references therein). Going to higher energies, we can study the V-A structure 
of the top quark couplings to the gauge bosons ( 7 , Z and W) [3]. The information about cou¬ 
plings can be obtained by using spin correlations: the top quarks are produced in certain spin 
states, as dictated by the top - 7 , 2’ couplings. Since the top decays before hadronization, 
the spin states of the top directly influence the angular distributions of its decay products. 
Simulations show that by analyzing kinematic variables of hnal state particles it is possible 
to measure the top-gauge boson couplings at the several percent level [ 1 ]. 

The subject of this paper is the evaluation of NLO QCD corrections to the top production 
and decay process above threshold at a Linear Collider. Previous work in this area includes 
the study of corrections to the top production subprocess, with virtual and soft gluons [4] as 
well as with hard gluon radiation ([5], [ 6 ], [7] are just some examples), and to the top decay 
subprocess with virtual and hard gluons together ([ 8 ], [9], [10]). 

Using these results, one can try to understand the top production and decay process 
by assuming that the intermediate top quarks are on-shell (narrow width approximation) 
and treating the subprocesses separately [11]. This assumption is usually reasonable; the 
result for the total cross section is valid up to terms of order Tt/rrit ~ 1%. However, the 
effects of hnite top width can be important in some differential cross sections; and if precision 
better than percent level is needed, the production and decay processes have to be considered 
together, by allowing the top quark to go off-shell. These corrections can be thought as being 
due to interference between production and decay, and are also known as nonfactorizable 
corrections. Some results concerning the nonfactorizable corrections for the top production 
process have already been published [ 12 ]; we present here a more complete analysis. 

We shall perform this computation and present the results at the parton level only. We 
assume that the issues related to jet reconstruction and identihcation can been solved, and 
the hnal state contains two W bosons ^ and two b quarks (the real gluon radiation case 
has been analyzed in [13], [14]). Even at this level, the complete computation of all the 
diagrams contributing to this hnal state (Born and next to leading order) is a very difficult 
task. Therefore, we shall employ the double pole approximation (DBA) which means taking 
into account only the diagrams which contain a top - antitop pair. 

Here, it is worth mentioning the strong similarities between the evaluation of QCD cor- 

^At the Monte Carlo level, we actually allow the on-shell IE’s to decay, either semileptonically or into a 
pair of massless quarks, but in the latter case we do not take into consideration QCD corrections to the W 
decay. 


2 



rections to the process of interest to us: 

e+e- ^bW+bW- (1) 

and the evaluation of QED corrections to the W pair production and decay process 
W~^ W~ —4/) at LEP. The issues which arise in the two computations are similar, because 
in both cases we are dealing with the production and decay of heavy unstable particles. Our 
treatment is largely similar to the one used for the electroweak process [15], [16], [18]. But 
there are some differences, both in the implementation of the DPA approximation and in 
the number and type of terms which contribute to the hnal result (the latter being due to 
the fact that in our case the intermediate off-shell particles are fermions, and not bosons). 
These differences will be pointed out in the course of our discussion. 

The outline of the paper is as follows. In sect. 2 we lay out in some detail the theoretical 
framework in which we perform our computation. This includes a short description of the 
DPA method, with examples of evaluation of NLO amplitudes in this approximation. The 
main point of this section is that the results for the amplitudes corresponding to interference 
diagrams are similar to results previously obtained for the W pair production process, while 
for the vertex corrections and fermion self-energy diagrams there are differences between the 
two cases. In order to facilitate comparisons with the on-shell approach, we also formulate 
our results in terms of corrections to the production and decay subprocesses and interference 
contributions. The gauge invariance of the total and partial amplitudes in DPA is manifest 
in this formulation. Sect. 3 contains some details on the design of the Portran code we 
used to obtain our numerical results, which are presented in sect. 5. A short discussion of 
a slighty different approach in evaluating the non-factorizable corrections (used in [12]) is 
presented in sect. 4. We end with the conclusions. 

2 Amplitudes in DPA 

Since it decays very fast, the top quark is not studied directly, but through its decay products. 
In top pair production at linear colliders, the relevant process is: 

e+e" ^ bW+bW- (2) 

At tree level, the diagrams contributing to this process can be split into 3 classes: diagrams 
which contain a top-antitop pair (Fig. 1), diagrams which contain either one top quark or a 
top antiquark (Fig. 2 and charge conjugates), and diagrams which do not contain any top 
(there are about 50 such diagrams). The evaluation of all these amplitudes can be performed 
using one of the automated tree level amplitude computation programs, like MADGRAPH 
([19], see also [20]). The computation of QCD corrections to all tree level diagrams, however, 
increases the degree of complexity by quite a lot, and is probably not feasible yet. 

Hence the need to use some methods which will simplify the computation, but will 
provide approximate results for the virtual corrections. One such method is the double pole 
approximation (DPA), which makes use of the fact that, looking for top pair production, 
we are interested in a specihc region of the hnal state phase space of (2). This region is 
dehned by the requirement that the invariant mass of the W b pair is close to the top mass: 
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Figure 2: Single top diagrams contributing to the process (2). 


pj, Pi ~ , where pt = pw+ + Pb and pi = pw- + Pi- In this region, the amplitudes 

corresponding to the top-antitop diagrams are enhanced by the two resonant propagators 
coming from the two intermediate top quarks"^: 


^ 9 — 9 9 —9 \^ / 

Pt - m pf - mi 

(the diagrams in (Fig. 1) are therefore called doubly resonant) and we can neglect the con¬ 
tributions coming from the singly resonant diagrams (Fig. 2) or nonresonant ones, which are 
reducedhj factors of (Tt/mt) or (Ft/rrit)^ with respect to the doubly resonant contributions. 

The real gluon corrections to the doubly resonant diagrams have been computed in [13]. 
The aim of this paper is the DPA evaluation of the virtual corrections to the top produc¬ 
tion and decay process. Some of Feynman diagrams contributing to these corrections are 
presented in Figure 3. These diagrams can be roughly divided into two classes : corrections 
to particular subprocesses - the vertex and fermion self-energy diagrams in Fig. 3 a) and 
b) respectively - and interference type corrections (Fig. 3 c) and d)). Strictly speaking, the 
vertex and self-energy diagrams also contribute to interference between subprocesses; but, 

"‘We use the notation m? = to? — imtTt 
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Figure 3: 


Feynman diagrams for virtual gluon corrections to top production and decay. 


for the sake of brevity, we shall refer to the diagram in Fig. 3 a) as the production vertex 
correction diagram, and so on. Also, in the following, we will denote the tree level amplitude 
(Fig. 1) by and the amplitude for the hrst order virtual corrections by 

= Mu + Mtb + Mu, + Mbt + Mu + M/,b ( 4 ) 

Here, the first three terms correspond to the three off-shell vertex corrections (which include 
in a suitable way the fermion self-energies, as described in section 2.4), and the last three 
terms come from the interference diagrams. 


2.1 Interference diagrams 

We start by discussing the evaluation of interference diagrams in Figs. 3c), 3d). Technically, 
this is one of the more involved computations we have to perform. The full evaluation of the 
t — b interference amplitude^: 


Mbi = u{b) 



d^k 1 ^ pb- Jk + mb , pt- ]k^mt 

k'^ + ie {pb — k)"^ — ml {pt — k)"^ — fhl 


X 


- h- 

{Pi + k) 


jk + mt 


— m 


— 7/^ 


- k + mt 
2 - 2 
pf - mi 


/w- v{b) 


( 5 ) 


is, for example, quite a difficult task. However, the only terms of interest to us in DPA are 
those which have resonances at the top and antitop quark propagator poles. This simplihes 
our task greatly. The doubly resonant terms can be extracted with the help of the following 
observation: if the virtual gluon in the loop is hard, then the quantity in brackets does 


®The Feynman gauge is used for the gluon propagator. 
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not have any singularity, and the overall resonant structure for this diagram is given only 
by the pole due to the antitop propagator; Af;,t(hard gluon) oc l/{pi — This means 
that any doubly resonant terms contribution to Aibi are entirely due to soft virtual gluons. 
Therefore, we can neglect the /c terms in the numerator of (5). Following [16], we shall call 
this approximation the extended soft gluon approximation (ESGA) 

With the help of the transformations: 

Y (A + = {- A + mb) A + ^ ‘Ipb (6) 


(- A + mt) 7^ = 7^ (A + mt) - "Ipip. -^Pif, 

(the term (A + mt) on the second line is neglected, since it would lead to a singly resonant 
contribution), the following result is obtained for the amplitude (5): 




Mbt{DPA + ESGA) = ^ x {-AbPt){pl-ml) 

47r 


1 


X 


d^k 1 _1_ 

k? + ie k"^ — 2kpb {pt — k)"^ — {pi + k)"^ 


ml 


(7) 


The result is proportional to the tree level amplitude ~ in the DPA, the virtual corrections 
due to interference factorize. The proportionality factor includes the scalar 4-point function 
(the integral on the second line of Eq. (7)) = V^{—pb, —pt,pt,0,mb,rht,mt)- ^ 

Let’s discuss shortly the resonant behavior of the DPA amplitude in (7). Apparently, 
the result for Aibi has a single pole at pj = mf (the other pole being canceled by the 
multiplicative term p'f —mf). However, if the top (or antitop) goes on-shell, the T>^ function 
acquires an infrared singularity (in the zero top width limit; this singularity is regularized by 
the top width). Since the infrared singular type terms have a logarithmic structure (this can 
also be reasoned from power counting arguments), this indicates that T>^ has the following 
behavior close to the top resonances: 


Vh ~ ailog(A - ml) + a2log(A 


ml) 


( 8 ) 


Here ai and 02 are terms which are hnite when either the top or antitop quark go on-shell. 
Formally, then, the overall resonant behavior of the interference amplitude AAbi in DPA is 

®In the standard soft gluon approximation, terms in the denominator of top quark propagators would 
also be neglected; we do not do this here for computational reasons (see also [16], [21]). 

^For the scalar one-loop integrals appearing here we use the following notation: 


I?°(P1,P2,P3, Wo,TOi,m2,TO3) 


1 

IVo Ni N2 N3 


£°{pi,P2,P3,PA,mo,mi,m2,m3,mA 
with the denominators : 


d^k 1 

No Ni N2 No Ni 


No = k"^ - ie, Ni = {k + pA - + ie, z = l,...,4 
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of type pole x logarithm: 


1 


Mbi^ Mo log(pf - ml) 


2 — 2 
Pt - 


rather than polexpole, as it is for the tree level amplitude or for the corrections to production 
or decay subprocesses. 

Using the same techniques, similar results are easily obtained for the other two interfer¬ 
ence diagrams. In the soft gluon approximation (and DPA): 

Ma{DPA + ESGA) = ^ M„ x (-4p,pj)(p? - in?) D,»s (9) 

Ma(DPA +ESGA) Mo x (-4psps)(pJ - - rnj) (10) 

where = V^{-pi, -pi, pt,0,mb,mt,mt) and 

^bb = -pt,PhPb,h,mb,Tht,mt,mb) ( 11 ) 


is the scalar 5-point function (here p, is the inhnitesimally small gluon mass needed for the 
regularization of infrared divergent behavior of 

We end this section with some comments on the numerical magnitude of interference 
terms. Since the resonant behavior of these terms is of pole x log type, it might be expected 
that they are less important numerically that the double pole terms. However, analytic 
expressions for the Vo function ([16], [17], [18]) show that, although the coefficients ai,a 2 
in (8) are finite when one of the top or antitop quark goes on shell, they will diverge when 
both particles go on-shell simultaneously: 

1 

cii{pl - ml) + C 2 i{pi - rhl) 

Therefore, the leading logarithms in the scalar 4 and 5-point functions will be enhanced by 
factors of order mt/Tf near the top, antitop quark mass resonances. 


2.2 Vertex corrections 


The results for the interference diagrams are completely analogous to results obtained in the 
W pair production computation. However, for the off-shell vertex and self-energy corrections 
diagrams, the results in the top case are different. Consider for example, the amplitude for 
the general vertex correction in Figure 4: 




CXs 

47r 


rd^k 1 


V 


^ 1 - Jk^mi 
{pi — k)"^ — ml 


r{Cv + Ca7^) 


- h- A + m2 
{P2 + k)"^ - ml 


Upon evaluation (and keeping only the vector part) the result can be written in terms of 
eight form factors, each of them multiplying a different tensor quantity: 

(Y 

= T^Cy Yi^F2 + ( a - mi)')^Fi -f 7^(- A - "^2) A 
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Figure 4: General vertex correction diagram. 


+ (j6i ^2 - m2)F8+p^Fi + -mi)p’(F8 + ...] (12) 

(expressions for the scalar form factors Fi,... ,F8 can be found in the Appendix). In the 
on-shell case, only the F 2 (electric dipole) and Fi (magnetic dipole momentum) form factors 
contribute. It might be expected that in the double pole approximation we can drop the 
other terms, too, since they have a zero at = mi (or ^2 = —^ 2 ) which will cancel one 
pole (or both) in the amplitude. However, the form factors themselves may have a resonant 
structure when the particles go on shell. 

Consider the top decay vertex correction. In this case, pi Pb,P 2 —^ —pt, and only 
four terms survive in Eq. (12); the corresponding form factors contain terms which are 
proportional to the scalar 3-point function: 


p ^ ro = /^ 1 1 _ 1 

J k‘^ {ph — ky — ml (pt — ky — ml 


f = 1,2,5,6 


which has a logarithmic resonant behavior: 


C[ 


tb 


\og{pl - my). 


Therefore, the contribution of f = 5, 6 terms to the top decay vertex correction is doubly 
resonant, although of type pole x log rather than double pole. Because these logarithms are 
not multiplied by large factors (as in the case of the interference diagrams), these terms can 
be expected to be numerically small; for consistency reasons it is still desirable to include 
them in the hnal result. 

Similar results are obtained for the correction to the antitop decay vertex (we keep the 
i = 1, 2, 3,4 terms in this case). In the case of the t — i vertex, though, both fermions are 
off-shell; as a consequence, there are no resonant logarithms when either the top or antitop 
quark goes on-shell, and we keep only the i = 1 , 2 terms. 

It follows that in the general expression (12) it is necessary to keep the terms which 
contain Fi to F^ (the F 7 and Fg terms can be dropped), and we don’t have factorization for 
the interference part anymore. This is different from what happens in the W pair production 
process, where the DBA factorization holds even for the vertex corrections. This difference 
between the two cases is due to the fact that in one process the intermediate particles are 
fermions, while in the other they are bosons. 
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Figure 5: Terms contributing to the renormalized vertex; the dots represent counterterm 
insertions. 


2.3 Renormalization and fermion self-energy 

For the renormalization of the ultraviolet divergences appearing in the evaluation of the ver¬ 
tex and fermion self-energy corrections we use the counterterm method. What this amounts 
to is replacing the bare vertex correction in Fig. 4 (which is UV divergent) with a hnite 
renormalized vertex correction (Fig. 5): 

= 6 T>‘ + mz,_ + h-it.A.Pi))^—+ —(-ffi 2 )(A) (13) 

2 Pi- mi 2 - p2-m2 

which also includes the contributions of the fermion self-energy diagram. 

The hrst two terms in Eq. (13) are what is usually dehned as the renormalized vertex. 
The last two terms are one half of the renormalized fermion and antifermion self-energy: 

S 2 (i 6 ) = S 2 (i>) - (± - m) 6 Z 2 - Am (14) 


where < 5^2 and Am are coefficients of the counterterms in the Lagrangian density (the ± 
sign corresponds to the particle/antiparticle case). In the above equation, S 2 stands for the 
bare fermion self-energy: 




f ^ /I X + m 

47r J {p — ky — m? 


(15) 


Upon evaluation, the result for Eq. (15) can be separated into a scalar and a spinorial 
component: 

S(j6) = (±/»-m)Sa(p2) -h mSfe(p2) 

With these notations, the renormalized vertex correction can be written as: 


= ff" + i ^zApA r" + i r- ^zArA 


with 


AZ2{p) = S„(p") + 


2 . ^ m'ZbijP') — Am 


A p — m 

The counterterm coefficient Am is hxed by the on-shell renormalization condition: 


(16) 

(17) 


f^ 2 i^ = m) = 0 ^ Am = mSfe(m^) 


(18) 
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Also, in the on-shell limit, 


AZ2ip)\^^m= Z.aip) + 


2 N , dJlbif i>) 


d jf) 


j/—^m j 


d i) 


= SZo 


(19) 




where for the last equality we have used the renormalization condition 

5^2 (i>) 


d i) 


= 0 




It is convenient to write the contribution of the fermion self-energy diagrams in a form 
similar to that of Eq. (12). Using the resummed top quark propagator: 




m] 


± ^ — m 

we obtain the following result for AZ 2 : 

AZ2{p)= f[S,(p2) + 2S,,(p2)] + 


p2 _ ^2 


( 20 ) 


S,..(p" 


m 


-(± ^ — m) 


( 21 ) 


with 


Sir(p^) = m 


Sfc(/) - Sfc(m^ 

p2 _ fj^2 


The term in square brackets in Eq. (21) will multiply the Born cross section. The term 
proportional to (± ^ — m) is identical to the like terms appearing in the expression for the 
vertex correction Eq. (12). Since is the part of the self-energy correction which would 
be infrared divergent on-shell (which means that it has a logarithmic resonant behavior 

Sir(/) ~ log(/ - m^) 
in the off-shell case), this term is also kept. 


2.4 Gauge invariance and corrections to particular subprocesses 

The partial amplitudes appearing in Eq. (4) can be directly related to Feynman diagrams 
and are straightforward to evaluate. However, as mentioned before, they cannot be directly 
identihed with corrections to particular subprocesses. For example, the top - antitop pro¬ 
duction vertex diagram (Figure 3a)) contributes to the correction to production vertex, as 
well as to interference between production and decay and even to interference between top 
decay and antitop decay, depending on when the top quark propagators are closer to the res¬ 
onances. Moreover, the amplitudes in Eq. (4) are not gauge invariant one by one, although 
their sum is. 

For purposes related to gauge invariance, and in order to be able to perform comparisons 
with the on-shell computation, it is desirable to decompose the total amplitude into gauge 
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invariant corrections to particular subprocesses, and interference between these. The aim is 
to rewrite Eq. (4) as: 




Mproii + Mtdpp + Midpc + + .^^h-fdec + ■^Zli-idpc 


( 22 ) 


with each term being gauge invariant by itself. 

To this end, it is necessary to decompose the amplitudes Mu, Mtb--- into parts which 
contribute solely to corrections to production, decay, or interference. This decomposition 
will be based on the top and antitop propagator structure of the matrix element. Following 
the prescription used in [13], products of propagators which go on-shell in different regions 
of the phase space can be decomposed as follows: 


1 1 _ 1 1 1 1 
0(pt) £>(p, - t) D„{p, - k) D{p,) D„(p, - k) D{pt - k) 


(23) 


with 

D{p) = — rfi^ , Dq{p — k) = {p — kY — p^ (24) 

In Eq. (23) the hrst term on the right hand side is considered as a contribution to the 
production process and the second one a contribution to the decay process. Furthermore, it 
is convenient to write the result in term of products of gauge invariant currents (in a manner 
similar to [16]). For example, the Mu amplitude can be written (using the extended soft 
gluon approximation): 


MuiESGA) 


47r 


Mf 


Ih r Ik) Dip,) 

iTT^ ' D{pi + k) D„{pp - k) Dip, - k) 


(25) 


where G^y{k) is the gluon propagator in an arbitrary gauge. By using the currents: 


3tdec 




D{pt 


Doipb - k) D{pt - k) 


(26) 


D,{p,+ k) 

the expression (25) can be written as: 


-‘^Pi D{Pt) 
DQ{pt + k) D{pt + k) 


MaiESGA) = Gppik) i poS - ifL" )],!: (27) 

where the hrst term in parentheses contributes to production-decay interference, and the 
second one contributes to decay-decay interference. 

We can similarly dehne the currents: 


D{Pi) 

Do{pi + k) D{pt + k) 


( 28 ) 


■ t, ^ ‘^Pt ‘^Pt DjPt) 

D„ip, - k) • D„(p, - k) D(p, - k) 
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and the amplitudes for the other two interference diagrams can be written like: 

M,i(ESGA) = G,Ak) ( - Ji" )*!„■' (29) 

M,dESGA) = ^mJ^ G^,(k) Jii" ■ (30) 

47r J 

Contributions to interference between subprocesses do not come solely from the mani¬ 
festly non-factorizable diagrams. The diagrams in which the gluon contributes to vertex or 
self-energy corrections (Fig. 3a)) also contain interference terms. Since the decomposition 
into purely vertex (or self-energy) corrections and interference corrections is not unique, we 
shall present our approach in some detail. 

The amplitude for the vertex correction diagram with off-shell particles can be written 
as: 


as fd'^k ^ A{pi,p2)+k^Ba{pi,p2)+ k^kdCap{pi,p-2) 

= S y 2- -((p. _ -7. (31) 


A corresponding on-shell approximation for this amplitude would be 

d^k _ ,,, A(p;,p^) + k^BM,P2) + k'^k/^Cc,f3{p'i,P2) 


mos = ^ 

An J in^ 


[—2p'ik + k^) {2p2k + kP‘) 


Iv 


(32) 


where and p '2 are some on-shell approximations for pi and p 2 {p\ = P 2 = m^)- Now, we 
can dehne the interference contribution through: 


(5r)i2 = (5r)- + (5r)r2^ 


intf 


(33) 


Note, however, that ((yr )°2 is not unique, since p[,P 2 are not unique; different choices for 
these momenta would yield different results for (5r)"|. The uncertainty which arises is, of 
course, of order p^ — rn^, so it can be neglected in the DPA. However, it allows us to choose 
the following dehnition for (^T)"!: 




An J in'^ _T>o(pi - k)Do{p 2 + k) 


r 


^k^Ba{pi,P2) + k°‘kdCal3{Pl,P2) 


7 


(34) 


D{pi - k)D{p2 + k) 

This choice means that the purely vertex correction (factorizable) part of the vertex diagram 
can be obtained by simply replacing the off-shell Co(pi,P 2 , 0 , mi, m 2 ) function appearing in 
the expression for (( 5 r)i 2 with the on-shell, infrared divergent function Co(pi,P 2 , P, \f^-i \f^)- 
Conversely, the interference part of the off-shell vertex correction diagram is: 




ds 

An 


d*k 

in"^ 


G,Ak) 


' M A{PI,P 2 ) 

^ D{p,-k)D{p2 + ky 


i 2 py r (- 2 p-) 

Do{pi - k)Do{p2 + k) 


(35) 
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with A{pi,p2) = (i^i + ^^2 + TTi2). For the tt production diagram, in DPA 


+ mt) Tti (- h + mt)Y ^ 

leading to 



Things are different for the decay vertex corrections, since we have doubly resonant 
contributions which are not proportional to the tree level matrix element. In the top decay 
case, the transformation: 

+ rhb) ^w+ (A + (M) ^w+ [2A + A + mt)] 

will lead to: 

= £ / 0 J.L" [(-vJjMo + M[ "] (37) 

where 

^ = n( \ m ^ A + hit) /w-yW] (38) 

D{pt) D{pi) L J 

In a similar manner, the interference term coming from the antitop vertex correction 

diagram can be written in terms of the currents Jidec i ~jprod^ the matrix element M(’ ^. 

Finally, the last diagrams to be split into on-shell and interference contribution are the 
top, antitop self-energy diagrams. Using the same approach as in the vertex case, we define: 


(az)A^ = 


CXg 

An 



t, h + mt ^ 
7 D{pt-k)^ 


^AV _11 

DoiPt -k)\) 


h+mt 

Dipt) 


^ fd^ 7,7 ‘^Pt 

Utt J tn^ ^ D,ipt-k)D,ipt-k) 


(39) 


where the quantity in the the curly brackets is the renormalized top self-energy, and the 
quantity on the second line is the on-shell limit of the quantity on the hrst line. This will 
lead to the following result for the interference contribution coming from the top self-energy 
diagram: 

Mu’ = ^ ^ [Ai* M, - M\- “] (40) 

and a similar one from the antitop self-energy diagram. 

Now we have all the pieces needed to write down the interference terms. The hnal result 
is: 

O' r ii*k r 

= ^ + M‘t - M[- “hlL 

+ jtdecJtdec-D^ o] (41) 
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It is easy in this formula to identify the production-decay or decay-decay interference terms. 
The currents: 


jprod jprod Jprod i jtdec Jtdec Jtdec ) jtdec Jtdec 


Jtdec 


(42) 


are conserved, and gauge invariant in DPA. Therefore, the total interference contribution as 
well as the interference between subprocesses parts are gauge invariant in the approximation 
used. 


3 Computational Approach 

Once a consistent scheme for evaluating the virtual corrections to the top production and 
decay process (1) has been set up (as described in the previous sections), the next step is 
the implementation of this scheme in a Monte Carlo program. In this section we give some 
details about the technical issues arising in the design of such a program, and how we choose 
to solve them. 

There are two types of quantities involved in the evaluation of the NLO amplitude: scalar 
quantities (form factors), which encode the contribution of loops, and spinorial quantities, 
built from Dirac spinors and operators. For example, the contribution coming from the ti 
vertex correction can be written: 

Mti = u{pb) + rrit) (5r(.g„(- ^i + m) /w-v{Pb) 

or, using the decomposition in Eq. (12) : 

= ? E [Ck Fl rr + Ca Ff- Tf 

4 = 1,2 ^ 

(for more details see [22]). 

Let’s start by discussing the evaluation of the scalar form factors F). Rather than compute 
analytic expressions for each form factor, we have chosen to evaluate them in terms of 
Passarino-Veltman (PV) functions [23]. This approach has the advantage that we have to 
compute only a few quantities which contain logarithms and dilogarithms: the Bq two-point 
and Cq three-point scalar functions (all the rest can be written as linear combinations of 
these functions). In turn, for evaluating the PV 2 and 3-point scalar functions, we use the 
FF routines [24]. 

To compute the amplitudes corresponding to the interference diagrams, we need to be 
able to evaluate the 4-point and 5-point scalar integrals in Eqs. (7), (9), (10). There are 
no published results (or routines) for the general (complex masses) 4-point scalar integrals. 
We have build such routines for the infrared hnite Vq function by using the general methods 
described in [25]. The results of these routines have been checked against analytical results 
in the soft gluon approximation published in [18]. 

The 5-point scalar function So has been computed by reduction to 4-point functions, 
following the recipe in [16]. The resulting infrared divergent 4-point functions have been 
evaluated using the analytic results published in [26]. 


(44) 


(43) 
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Some comments on the treatment of the top width are needed here. One way of evaluating 
the scalar form factors in Eq. (44) is to compute the gluon integrals in the zero top width 
limit and introduce the hnite width only in terms which are divergent on-shell (that is, 
replace mf with = mf — imtVt in terms like log(p^ — rnf)] see for example [16]). The 
difference between this result and the one obtained by using the complex top mass in all 
the terms is of order Tt/mt, therefore at about 1% level. This would be acceptable if the 
radiative corrections would be small with respect with the tree level result (as is the case for 
the W production process), but in our case it turns out that the one-loop QCD corrections 
are of the same order of magnitude as the tree level result®. Therefore, order % terms are 
important. Since in the case of real gluon radiation the top width appears in all terms, for 
reasons of consistency we need to keep the width in all terms in the evaluation of the virtual 
corrections too. 

The other elements needed in the evaluation of the amplitude (44) are the spinor sand¬ 
wiches Tj. We compute these quantities using spinor techniques, as for the real gluon radi¬ 
ation case. Since this part of the computation is quite complex, and hence prone to errors, 
we have two different ways of performing it. In one approach, we express the Tj’s in terms 
of basic spinor products u{pi, Si)u{pj, Sj)] this is the more involved case (in terms of the 
work done by the programmer), geared for implementation in a Fortran routine, and which 
allows fast computation. The other approach uses C-I--I- routines which allow the automated 
evaluation of general spinor sandwiches like 

h(p, s)(A +mi )(/>2 + 1712 ). ..u{p',s) 

(To this purpose, we have constructed classes that describe < bra\ and \ket > spinors, and 
operators of type /f)i ± mp in turn, these classes use the basic classes - 4-vector, complex 
number - dehned in the Pandora event generator [27]). This method allows easy evaluation 
of Tj expressions (again from the programmer’s viewpoint) and is much more resistant to 
programming errors; but the computation is slower than in the previous method. Therefore, 
the main use of the results obtained from the C-I--I- routines is to check the Fortran results. 

4 On-shell DPA 

The issue of interference effects in the production and decay of heavy unstable particles has 
been the subject of extensive studies in the past decade. One of the main results is the 
theorem, due to Fadin, Khoze and Martin [28], which states that these interference effects 
are suppressed (see also a more recent discussion in [29]). A stronger version of this theorem 
[30] claims that NLO interference effects cancel in inclusive quantities up to terms of order 
aV/M. As in [30], it is possible to dehne and use a framework for the computation of 
interference corrections in which the total interference contribution to inclusive quantities is 
zero. 

In order to facilitate a comparision with our results (where the interference contributions 
do not cancel completely), we shall discuss this alternative approach (which we shall call on- 

®The reason the total QCD corrections are of order 10 - 20% is because of large cancellations between 
the virtual corrections and soft gluon real corrections. 
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shell DPA) in some detail. Results obtained using on-shell DPA have been obtained for the 
W pair production case at LEP II ([16], [18]), as well as, more recently, for the evaluation of 
interference (non-factorizable) corrections to the top production and decay process at 
and hadron colliders [12], 

The relevant features of this approach are two: first, the amplitudes for corrections to 
subprocesses are computed in the on-shell approximation for the top quarks. For example, 
the correction to the production process can be written in terms of the on-shell amplitude: 

= E ^ “■(»)) -^aI* ^ »'■+) My(t ^ m-) (45) 

A,A' 

where A, A' are the spins of the top quarks. The difference between the above amplitudes and 
the ones used in this paper (Eqs. (22), (34)) is due to non-doubly resonant terms, therefore 
acceptable in DPA. 

The other characteristic feature of the on-shell DPA method is that the interference due 
to real gluon radiation is computed by using a semianalytic approach. This approach rests 
on the observation that interference is due mainly to gluons of energies of order P^; therefore, 
the (extended) soft gluon approximation can be employed for the evaluation of interference 
terms. 

There are two stages where this approximation comes into play. First, it is used at the 
matrix element evaluation level. For example, consider interference between the diagrams 
where the gluon is radiated from the bottom quark and from the antitop quark: 


d(TlhiPb,Pw, ■ ■ ■ ,k) ~2Re 


E -Mr (Mir 


(46) 


= |Afo(p&,Pw, • • OP 2Re 


^PtPb Pt - ml _ 1 _ 

2kpi, + ie {pt + ky — frit {Pt + ky — fht* 


The second stage is the treatment of the hnal state phase space. In the soft gluon 
approximation, it factorizes: dQb,w,...,g = dQb,w,... x dQg , and the integration over the gluon 
momenta is performed separately: 


d(T^ibiPlPw^---) 


l-Mo(Pfe,PvK, • • •)p 


X 



d^k ^ptPb Pt~^t 1 

27ra; 2kpb + ie {pt + ky — fht {pi + ky — fht* 


(47) 


where p'b^p'w) ■ ■ ■ given by a suitable projection of the off-shell momenta Pb,PW: ■ ■ ■ onto 
the on-shell phase space (for an example of how this projection might be accomplished see 
[31]). In the above equation pt = p], -f p(y+, Pt = p'l + and the integral over gluon 
momenta is allowed to go to inhnity (since hard gluons contribute nonresonant terms to the 
result). 

The quantity on the second line of Eq. (47) can be evaluated analytically, through 
methods similar to those used to evaluate the virtual 4-point functions. We will not give the 
results here (they can be found in [16], [18]), but there is an important comment to make. 
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Using this procedure to compute the real gluon interference, the total interference obtained 
by adding the virtual diagram contribution (Eq. (7)) to the above result and integrating 
over the top invariant mass parameter is zero [30]. This cancellation works also for the 
other interference diagrams; therefore, in this approach, the contribution of non-factorizable 
corrections is zero to the total cross section. However, this result depends on two things. 
First, it requires an inclusive treatment of real gluon radiation, with phase space integration 
extending to inhnity. Second, both the virtual and the real interference terms have to be 
treated in the soft gluon approximation. 

But, is the use of the soft gluon approximation justihed in this case? At the amplitude 
level (Eq. (46)), the answer is yes; the relevant gluon energy, being of order T^, is much 
smaller than the other momenta involved. However, this approximation does not seem to 
be acceptable for the phase space factorization stage of the above approach. Here, problems 
might arise when projecting the off-shell momenta onto the on-shell phase space. The reason 
is that there is no single way to perform this projection; therefore, in the determination of 
the on-shell momenta ... there is an uncertainty of the order of the gluon energy, or 

Tt. Now, being close to the top resonances, we are in a region of the phase space where 
the cross section varies greatly over a range of energy of order Tt (due to the top quark 
propagators); therefore such an uncertainty is not acceptable. 

To illustrate the dependence of the result for real gluon interference on the choice of the 
on-shell momenta pj,, p(y,..., let’s presume that instead of projecting pf, into we also take 
into account the gluon momentum: ph + k ^ p[ (physically, this might be justihed by the 
inclusion of the gluon jet in the bottom quark jet). Then, Eq. (47) becomes: 


dcr^biPlPw:---) 


Q;, 


\Mo{plp'w,---)r -Re 

71 


d^k AptPb _ 1 _ 

27ra; 2 kpb + ie {pi + kY — rfY* 


(48) 


The result for the above expression is different from the result for Eq. (47), and the difference 
contains doubly resonant terms. Therefore, in the on-shell DBA approach, the result for 
the interference terms depends on the choice of the implementation of the phase space 
factorization. A discussion of this dependence for the W pair production case can be found 
in [16]. 


5 Results for virtual corrections and the total cross 
section 

In this section, we present some results on the total cross section for the top production 
and decay process at linear colliders. We take into account the virtual corrections as well 
as contributions coming from real gluon radiation. Furthermore, we study the effect of 
interference (nonfactorizable) terms on invariant top mass distributions and we perform 
comparisons with results previously published [12]. 

In obtaining the results presented in this section, the following set of parameters is used: 

mt = 175 GeV, = 0.1, r° = 1.55 GeV, T^ = 1.42 GeV, 
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where r° is the top width at the tree level, while includes the NLO QCD radiative 
corrections. 

We start by looking at the total cross section for our process. Table 1 presents results 
for the following quantities: 

• (To : cross section for the tree level process (2); computed in the on-shell (narrow width) 
approximation, using the zero-order top width F^. 

• (T°^ : cross section for the NLO process in the on-shell approximation; computed using 
the NLO top width F^. 

• : the main (factorizable) part of the DPA approximation to the NLO process. 
This quantity contains corrections to production and decay as dehned in section 2.4. 

• : the interference (non-factorizable) part of the DPA approximation to the NLO 
process, as dehned in section 2.4. 



360 GeV 500 GeV 1000 GeV 

ao{pb) 
(t“ {pb) 
ai^'ipb) 

0.386 

0.700(1) 

0.676(2) 

-0.032(1) 

0.570 

0.660(2) 

0.664(2) 

-0.0116(3) 

0.172 

0.1839(7) 

0.1920(8) 

-0.0061(2) 


Table 1: Total cross sections for top production at linear colliders with no cuts on phase 
space (numbers in parantheses are errors due to numerical integration). 

We present results for three values of collision center-of-mass energies: 360 GeV, just 
above the tt production threshold, 500 GeV, the most common value used in linear collider 
studies, and 1 TeV, which is relevant for higher energy machines. The 360 GeV result is 
probably not good; this close to the threshold, resummation of large logarithms appearing 
in the tt interaction is needed [2]. However, it is interesting to see the magnitude of the 
nonfactorizable corrections at hxed order in this energy range. 

The NLO cross sections contain contributions from the virtual corrections as well as 
from real gluon radiation (evaluated as in [13]). The phase space splicing method is used 
for the treatment of infrared singularities. The value of the technical cut which separates 
the infrared from the real gluons is e = 0.1 GeV; the results are independent of the choice 
of this parameter®. No physical cuts have been imposed on the hnal phase space. 

The results in Table 1 prompt several comments. First, the difference between a°^ and 
(t{“‘^* is due to non-doubly-resonant terms, therefore it could be expected to be small. This 
is indeed the case at 500 GeV; but at 1 TeV, this difference is about 4% of the cross section. 
The reason is that in obtaining these results, we have integrated over the complete kinematic 

^Taking its value much smaller than the top width allows the use of the soft gluon approximation and 
phase space factorization in the evaluation of the infrared integrals. 
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range available for the top quark invariant mass (that is, mb+mw < yPt < W—{mb+mw)), 
so the cross-sections include contributions from regions of the phase space where the top 
quarks are far off-shell and non-resonant terms are important. 



360 GeV 500 GeV 1000 GeV 

(pb) 

a^r^pb) 

0.682(1) 

0.670(2) 

-0.034(1) 

0.627(2) 

0.629(2) 

-0.0068(2) 

0.1742(7) 

0.1781(5) 

-0.0017(1) 


Table 2: Total cross sections for top production at linear colliders with cuts on the top, 
antitop invariant mass. 


In Table 2 we present the cross section results obtained with a cut on the t, t invariant 
mass IY^, — mt I < 15 GeV. The difference between the two results for the main terms 

and is smaller in this case. Note that, since in the on-shell approach 
in Table 1 and 2 contains a factor which simulates the effect of cuts (either from kinematic 
constraints or imposed ones) on the top invariant mass. Note also that these cuts are not 
imposed ad hoc, but they arise rather naturally in the process of dehning a t,i production 
event; it makes sense to require that the reconstructed mass of the b, W pairs is close to 
the top mass in the dehnition of such an event. In this context, it is also worth noting that 
the contribution coming from the phase space region where either the t or t is far off-shell 
(more that ten times the width) is quite sizable (around 5% of the total cross section for 
center-of-mass energies greater than 500 GeV). 


Another quantity of interest is the differential interference cross section as a function of 
the top invariant mass. Even if the total interference contribution to the cross section is 
small (at about 1% level), it can have larger effects in differential distributions since it can 
be positive in certain regions of the phase space and negative in others. In particular, it can 
be important in the reconstruction of the top invariant mass; since tends to decrease 


as ypl increases, the net effect would be to shift the position of the Breit-Wigner peak to 
smaller invariant mass values. This effect can be quantihed by the following equation: the 
shift in the mass is 


AMt 


d5nf \ 

dMb) 


r? 


Mt=mt 


(49) 


where Mt = ypf, and 6 nf is the ratio of the non-factorizable (interference) part of the cross 
section to the Born cross section: 

^ daT'^ 

“ dao 

In Figure 6 we present the differential distribution for the relative non-factorizable cor¬ 
rection Snf{Mt) at center of mass energy 500 GeV. The dashed line is the result which takes 
into account the full interference corrections in Eq. (41); the solid line is obtained by taking 
into account only the terms proportional to the Born amplitude. Note that, although the 
contribution of the Mi terms in Eq. (41) to the total cross section is very close to zero. 
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nf 


0 


- 0.02 


- 0.04 



170 175 


180 Mt (GeV) 


Figure 6: The relative nonfactorizable correction to the invariant mass distribution; the solid 
line is the contribution of terms proportional to the tree level amplitude, while the dashed 
line contains also the Mi terms in Eq. (41). 


they have a sizable effect on the differential distribution in Fig. 6. Using Eq. (49), we 
conclude that the shift in the position of the peak in the top invariant mass distribution due 
to interference effects is very small (of order of a few MeV - in agreement with [12]). 

The results in Table 1 and 2 indicate that the contribution of interference terms to the 
total cross section is of order 1%, in agreement with the Tt/rrit order of magnitude expected 
from naive arguments. However, it is not zero, as implied by results presented in [12], which 
use the on-shell DPA method. We have argued in section 4 that this difference is due to the 
way in which the radiation of real gluon with energies of order F^ is treated. In Figure 7, we 



Figure 7: Real gluon interference: the abt term. The solid line corresponds to the semianalytic 
approach; the dashed line is obtained through numerical evaluation with Mt = \Jplw'-, fhe 
dotted line is obtained through numerical evaluation with Mt given by Eq. (50). 
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present the results for the real interference between the diagram where the gluon is radiated 
from the bottom quark and the diagram where the gluon is radiated from the antitop quark. 
The solid line is the result of the semianalytical approach described in section 4. The other 
two lines are the result of the exact off-shell computation (where the integration over the 
gluon momenta is performed numerically). The two lines correspond to two different ways in 
which the gluon momentum is treated in the reconstruction of the invariant top mass. For 
the dashed line, the gluon momentum is ignored in the top mass reconstruction: = \J'Pbw 

Note that in this case, the result is quite close to that of the semianalytical computation, 
which is natural, since the gluon momentum is treated in both cases the same way. 

To obtain the dotted line, we have followed a more realistic approach, in which the gluon 
is included in the top mass reconstruction if it happens to be radiated close enough to the 
top quark: 


Mt 


\jplwg if cos6'tg < 7r/3 
\JpIw otherwise 


(50) 


Although the total cross section is the same as for the other exact evaluation case, the 
differential cross section differs by quite a bit. 

The total cross section corresponding to the interference term presented in Figure 7 has 
the value = 0.121 pb for the semianalytical result, and ati = 0.124 pb for the numerical 
one. Note that, since the diagram set contributing to this interference term is not gauge 
invariant, this result has no physical meaning by itself. However, from these numbers we 
can get some insight concerning the evaluation of interference corrections. First, note that 
the contribution of this single diagram is much bigger (about two orders of magnitude) than 
the total result for the interference terms. This means that there are large cancellations 
taking place between the real and virtual interference contributions. This is quite natural, 
in accordance with the discussion in section 4; however, this also means that small uncer¬ 
tainty (order percent) in evaluating one of this contributions (the one coming from the real 
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Figure 8: The relative nonfactorizable correction to the invariant mass distribution; com¬ 
parison between the semianalytical (dashed line) and the numerical (solid line) approach. 
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gluon interference, for example) can lead to large uncertainties in the evaluation of the total 
interference contribution. 

In Figure 8 we present the comparison between the relative non-factorizable corrections 
computed in the semianalytical approximation (dashed line) and the complete off-shell ap¬ 
proach (solid line). For the purpose of this comparison, we consider only the terms propor¬ 
tional to the Born amplitude in the exact computation, since only these terms are taken into 
account in the semianalytical approach. Note that the total interference cross section inte¬ 
grates to zero in the latter case, and, as discussed above, the complete off-shell distribution 
contains contributions that do not cancel in the total cross section. 


6 Conclusions 

In this paper, we have discussed in some detail the evaluation of next to leading order 
QCD corrections to the top production and decay process at a linear collider. Since a full 
computation of the NLO amplitudes contributing to the process e+e“ —b b W~ is not 
feasible, the double pole approximation (DPA) has been employed. In our case this means 
taking into acconnt only the diagrams which contain two intermediate top qnarks. Unlike 
most of the previous treatments, we allow for the two top qnarks to be off-shell, and include 
also the corrections dne to interference between the top prodnction and decay processes. 

A parallel is drawn between our computation and the DPA evalnation of QED corrections 
to the W pair prodnction and decay process at LEP. While there are many similarities 
between the two compntations, there are also some differences; the most striking one is 
that there are nonfactorizable (interference) corrections no longer proportional to the Born 
amplitnde in the top qnark case. Althongh the contribntion of these terms to the total cross 
section (or other inclusive quantities) is consistent with zero, they contribnte to differential 
distribntions. 

Previous results obtained using a semianalytical approach state that the interference 
terms cancel ont completely in inclusive quantities (like the total cross section). We discnss 
the evalnation of the real glnon interference terms nsing analytic methods and we point ont 
the shortcomings of this approach. The total magnitnde of nonfactorizable corrections in 
our computation is fonnd to be of order 1% of the cross section. We also present resnlts 
for the total cross section and the relative nonfactorizable correction to the top invariant 
mass distribntion. The effect of nonfactorizable corrections on the top mass reconstrnction 
is fonnd to be very small. 

The results presented in this paper are relevant at collision energies above the top - 
antitop prodnction threshold. Since we are interested in differential distribntions of hnal 
state kinematic variables, the approach used to obtain cross section distribntions is that 
of nnmerical simulations. The amplitndes are evalnated nsing spinor techniqnes, and the 
integration over the hnal state variables is performed nsing Monte Carlo techniqnes. This 
approach has the added advantage that it allows for the inclusion of experimentally relevant 
selection criteria on the hnal state phase space. 

The calculation presented here is entirely at the parton level. For more realistic sim- 
nlations, it is necessary to take into acconnt initial state related issnes, like initial state 
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radiation (ISR), beam energy spread and beamstrahlung. The hadronization of the hnal 
state partons also has to be modeled. Further extensions of interest include taking into 
account singly-resonant diagrams (tree level as well as NLO in the on-shell approximation) 
and the evaluation of electroweak corrections. 
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Appendix 


In this appendix we shall give some details about the evaluation of vertex corrections 
and fermion self energies. For a more extensive description, we refer the reader to [22]. 

The amplitude for the general vertex correction in Figure 4 can be written as: 


(5F^ 


47r 


rd^k 1 


V 


/>!- jk^mx 
{pi — ky — fhi 


7^(C'y + CaI^) 


- h- A + m2 
{P2 + ky - ml 


(51) 


Upon evaluation of the integral, the result can be written in terms of two sets of form-factors; 
one for the vectorial part of the vertex correction, one for the axial part. The number of 
form factors needed depends on the specihc constraints on the process; in the general case, 
when the momenta pi,P 2 are off-shell, eightx 2 form factors are involved: 


■’r" = ^ E (CvFTtP + CAFf'T,''") (52) 

The dehnition of these form factors depends on the choice of the spinorial elements in 
terms of which the result is written. We dehne: 

5F0 = X (53) 



[ 

Pi 



fY 

+ 






FY 

+ 

(A 

- mi) 

Pi 



FY 

+ 

(A 

- mi) 




FY 

+ 



Pi 

(- A 

- m2) 

FY 

+ 




(- A 

- m2) 

FY 

+ 

(A 

- mi) 

Pi 

(- A 

- m2) 

FY 

+ 

(A 

- mi) 


(- A 

- m2) 

FY 

] 


for the vectorial part of the vertex correction; for the axial one, replace 7 ^ with 7 ^ 7 ® 
in the expression above. This definition has the advantage that when the particle i is on- 
shell, the terms which contain the ± drop out. Also, we have made use of the fact 
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that, if (5r^ multiplies in the full matrix element can be thought of as the polarization 
vector of the gauge boson in diagram 4), then (pi + P 2 )^A^ = 0. 

We shall evaluate and write the results for the form factors in terms of Passarino - 
Veltman functions (for the dehnition of these see below): 

= 4[ mi{Ci2 — Cii + 6*23 ~ C 21 ) — rn2{Ci2 + 6 * 23 ) ] (54) 

= -2[ 2pip2(C'o + Fn) + 2(^24 - 1/4) + 2Bl^ - 1 

—mim2Cii + Pi(—Cii + C 12 ) — P2^i2 ] 

=- 4(^0 + 2C'n - ^12 + ^21 - ^ 23 ) , F/ = 27712(^0+ C'n) 

F ^ = — 4 ( 6*0 + Cii + 6*23 + ^*12) ; Fq = 27771(6*0 + 6*11) 

= 0 , F^ = 2(6*0 + C-n) 

In the computation of the axial form factors, 7 ® can be shifted to the right in Eq. (51). 
In the result for the form factors, this amounts to changing the sign of m 2 , m 2 —m 2 , and 
multiplying the F 5 ,..., Fg with (—1). Then : 

Fi^ = F/" + 87772(612 + 623 ) , F^ = - Amim 2 Cii (55) 



The result (21) for the fermion self-energy corrections can be similarly written in terms 
of Passarino-Veltman functions: 

Sa(/) = ^(l + 2Fl(p^r772)) (56) 

4:71 

S-(p') = [4AFo(p^ 7772) +4AFl(p^ 7772)1 

47r — m^ L J 

with 

ABnip'^, rfi^) = Bnip^, m^) - Bn{rfi^, m^) , 77 = 0,1 
If we further dehne the Aq, Ai form factors through: 

AZ 2 {p) = 2 ^ [Ao(p2) + Ai(p2)(± ^ - m)] (57) 

471 L J 

the renormalized vertex correction in Eq. (16) can be obtained by making the following 
redehnitions of form-factors in Eqs. (54), (55): 

F^'^ ^ + Xo{pl) + Xo{pl) (58) 

F^F ^ F^^^ + X,{pI) 
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These general results are easily translated for the specihc cases which appear in our 
computation. For corections to the production, top decay and antitop decay vertices we 
have, respectively: 

Kit = = Pu P2 = Pi, mi = 7712 = m) (59) 

fK = K^{pi = Pb, P 2 = -pt, mi = mb, m 2 = mt) 

Kit = ^KiPi = -Pi, P2 = Pb, mi = mt, m 2 = mb) 


The statndard dehnition for the (3-point) Passarino-Veltman functions is: 


^ f<K_ __ 

J (/c^ — mf) {{k + pi )2 — m2) {{k + pi + ^2)^ — m\) 

with the scalar functions Cij{p‘l,p 2 ,m‘l,ml,ml) given by: 

= p^iCii+p^2Cl2 

= p>ip’iC2l+P^2P2C22 + {p^lP'^2+P^2Pl)C23+g^''C2, 


(60) 


(61) 


The expressions for the reduction of Cij functions in terms of the scalar one-, two- and 
three-point integrals A^, are standard, and we do not give them here (see [22] or [32]). 
The C function used in Eqs. (54), (55) are dehned by: 


J {kP'+ ie) {{k — pi)‘^ — mi) {{k + p2)'^ — m2) 


(62) 


= -p‘iCii + {p2-pi)^Ci2 (63) 

= P 1 P 1 C 21 + {P2 - Pl)^{P2 - Pl)‘'C22 - [PliP2 - Pi)"' + {P2 - Pl)>l] <^23 
+ (7^"C'24 

therefore: 

Co,Cij = C^,Cij{p\,{p2-pi)‘^,Q,ml,ml). (64) 

Moreover, the function in Eq. (54) is given by: 

Bl^ = B\{pi+p2f,ml,ml) (65) 

and in Eqs. (56): 

Bn{p^, m^) = Bn{p^, 0, m^) , forn = 0,1. (66) 
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